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The interfacial coupling of two materials with different ordered phases, such as a superconductor 
(S) and a ferromagnet (F) is driving new fundamental physics and innovative applications. For 
example, the creation of spin-filter Josephson junctions and the demonstration of triplet supercur¬ 
rents have suggested the potential of a dissipationless version of spintronics based on unconventional 
superconductivity. Here we demonstrate evidence for active quantum applications of S-F-S junc¬ 
tions, through the observation of macroscopic quantum tunneling in Josephson junctions with GdN 
ferromagnetic insulator barriers. We prove a clear transition from thermal to quantum regime at a 
crossover temperature of about 100 mK at zero magnetic held injunctions which demonstrate a clear 
signature of unconventional superconductivity. Following previous demonstration of passive S-F-S 
phase shifters in a phase qubit, our result paves the way to the active use of spin hlter Josephson 
systems in quantum hybrid circuits. 

Introduction 

Superconductor (S)/Ferromagnet (F) interfaces are a special class of hybrid systems where different ordered phases 
meet and generate new nanoscale phenomena and new forms of global order [TJ [5]. A ferromagnetic barrier in a 
Josephson junction (JJ) generates novel physics [3l |4], and represents a technological key for advances in weak super¬ 
conductivity, spintronics and quantum computation[S]-[7]. Recent interest in ultra-low-power, high-density cryogenic 
memories has spurred new efforts to simultaneously exploit superconducting and magnetic properties so as to create 
novel switching elements having these two competing orders[8]. S-F-S junctions are expected to shed light on several 
aspects of unconventional superconductivity, including transport through spin aligned triplet Cooper pairs [H[SHU]. 
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The NbN-GdN-NbN junctions investigated in these experiments are spin filter devices[TS] with an unconventional 
predominant second harmonic current-phase relation (CPR)[T6]. By changing the thickness of the GdN ferromagnetic 
insulator (FI) barrier it is possible to change its magnetic properties and hence the spin hlter efficiency fSFEl [T51ll7j . 
In this paper we report measurements on NbN-GdN-NbN junctions for very low and very high SFE values, from 
almost a few percent up to 90%, with an intermediate value of 60%. By increasing the thickness of GdN we also 
obtain junctions with lower values of the IcRn product, where Ic is the critical current and Rn the normal state 
resistance. Following the previous demonstration by Feofanov et al.[T8] of passive S-F-S phase shifters in a phase 
qubit, the experiment reported here demonstrates the active quantum potential of S-F-S JJs via the occurrence of 
macroscopic quantum tunneling (MQT) in spin filter devices. 
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Results 

Transport properties of spin filter Josephson junctions. 

The hallmark of the spin-filter effect is the decrease in resistance in the R vs T curves below the ferromagnetic 
transition, as shown in Fig. la {Tcurie — 33 K for GdN). SFE is defined as the percentage difference in the tunnelling 
probability for up/down spin electrons owing to the difference in barrier heights of the up/down spin channel in the 
FI caused by exchange splitting, so that 100% SFE corresponds to pure tunnelling of one spin sign. See Methods 
for details concerning the calculation of SEE. Table 1 collects the parameters of the measured junctions. All the 
NbN-GdN-NbN JJs present hysteresis larger than 90% in the current-voltage (IV) characteristics, in a wide range of 
Josephson coupling energies Ej = (where Ico is the critical current in absence of thermal fluctuations and 

4>o = /i/2e is the quantum flux). Eor all the measured devices the junction area is about 7^m x 7/im ~ 50 /im^. 

In Fig. Id, IV curves, measured as a function of temperature T for the junction with the highest SFE, are reported. 
The dependence of Ic as a function of 77 at T = 4 K is reported in Eig. Ic. The blue curve is the first measurement of 
the magnetic pattern, after nominal zero-field cooling. Ic{H) has then been measured both in the downward direction 
of the magnetic field sweep and in the upward direction (black and red curves, respectively). The black and red curves 
show a distinctive shift of the absolute maximum of Ic from — I mT to -I-1 mT respectively arising from the hysteretic 
reversal of the FI barrier[T^. The period of Ic{H) in non-spin filter junctions for the same geometry is twice as large 
(3 mT), pointing to a largely predominant second harmonic in the CPR in spin filter JJs, as discussed in Ref.|I 6 j. 

Measurements of switching current distributions. 

We have studied the escape rate of the superconducting phase (p as a function of temperature and for different values 
of the magnetic field, through standard switching current distribution (SGD) measurements. SCDs have been 
performed for all the samples reported in Table ||| According to the Resistively and Gapacitively Shunted Junction 
(RCSJ) model [231 US], for a JJ with a conventional CPR the dynamics is equivalent to that of a particle of mass 
moving in a washboard potential U{(p) = —Ej{cos(p + 'yp). The particle mass is given by = C'((()o/27r)^ with C 
the capacitance of the junction. The normalized bias current 7 =//Ico determines the tilt of the potential. 

The second harmonic component in the CPR, / = /isini^ -|- / 2 sin 2 (/j, leads to a modified washboard potential 
U{p) = —Ai(cos (/3 -I- I cos 2(/3 -|- El =hli/2e, which may assume the form of a double well for values of 

g=l 2 /Ii larger than 0.5 |26lE5] (see Fig. 2a). The presence of two wells in the washboard potential may result in 
the observation of two critical currents in the IV characteristics, since when tilting back the washboard potential the 
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phase particle may be retrapped in one of the potential wells with finite probability |27j: the case of JJs” with g 
< —1/2 has been recently studied [28l I29j . Measurements of two-well distinguished critical currents constitute a very 
direct criterion to estimate the g factor [?n [25] , Although spin filter JJs have a strong second harmonic component 
in the CPR [Ibj we have not found evidence of two critical currents, and hence the case of (p junction with negative 
values of g is not considered in this work. 

The sin2(p term in the CPR on average lowers the barrier height of the washboard potential without significantly 
altering the asymptotic expression of the potential barrier for 7 close to one. The height of the potential barrier 


IS 


given by: AU{j,g) = ^Eja{g) {1 —, where a{g)‘^ = plasma frequency is given by: 


9 ) = ^ 2 /a( 5 ) 2 i/ 2 e/co/SC'(l - 7 ^ In addition, the second harmonic component does not modify the power 
law of AU. As a consequence in the thermal regime the standard deviation a of the SCDs is expected to scale as 
as in the case of the standard CPR |20[ [^ . 

The motion of the particle is subject to damping given by where Q = ujpRC is the quality factor and R is the 
resistance of the junction. When the bias current is ramped from 7 = 0 to 7 < 1, the junction is in the zero voltage state 
in absence of thermal and quantum fluctuations, corresponding to the particle being localized in a potential well. At 
finite temperature the junction may switch into the finite voltage state for a bias current 7 < 1. This corresponds to the 
particle escaping from the well either by thermally activated processes or by tunneling through the barrier potential 
(see Fig. 2b). In the thermal activation (TA) regime, the escape rate for weak to moderate damping {Q > 1) is 
determined by[30] P* = exp where the thermal prefactor is Ot = 4 (1 -|- QkBT/1.8AU)^^^ + 1 [3T] . 

The escape rate will be dominated by MQT at low enough temperatures [35]. For Q > 1 and 7 close to 1 the escape 
rate in the quantum regime is: F^ = Og^exp where Og = (8647rA[//Iia;p)^^^. The crossover temperature 

between the thermal and quantum regimes is given by [33] 
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The experimental probability density of switching is related to the escape rate through the following equation]! 




f' nn 
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dR 
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where AI/At is the current ramp rate. 
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The measurements have been performed in a dilution refrigerator with reaches a base temperature of 20 mK. A full 
description of the apparatus is given in detail in the Methods section. The bias current of the junction is ramped 
at a constant sweep rate AI/At = 2 mA s“^ and at least 10^ switching events have been recorded using a standard 
technique [3S] . 

Fig. 3a shows a set of SCDs as a function of temperature for the high SFE JJ reported in Fig. 1. The thermal 
behavior of the SCDs is typical of underdamped JJs and the standard deviation a, which is proportional to the width 
of the switching histograms, increasing with temperature as expected. Fig. 3b shows the SCDs measured below 
IK (black circles). The dependence of the standard deviation a on temperature is reported in Fig. 4a (right axis), 
along with the thermal behavior of the mean value of the SCDs, Imean, below 0.5 K (left axis). When decreasing the 
temperature, Imean increases while a decreases and both saturate at a crossover temperature of about 100 mK. Below 
this crossover the histograms overlap and the escape process is no longer regulated by thermal fluctuations, indicating 
the transition to the MQT regime [201 [3T]. 

A further confirmation of the observation of MQT comes from measurements of SCDs in magnetic field. The 
behavior of o'(T) aA, H = 1.1 mT is shown in Fig. 4b. At this value of the magnetic field, which reduces Ic to half of 
the value measured at zero held (see the blue squares in Fig. Ic), lower values of a have been measured and Tcross 
is reduced by a factor down to about 70 mK, in agreement with MQT theory [55]. In both cases of OmT and 
1.1 mT, Tcross has been determined by the intersection of the T'^/^ curve in the TA regime (dashed green lines in 
Fig. 4b) and the mean value of a in the MQT regime (black full lines in Fig. 4b). The measurements in presence of 
magnetic held prove that the battening of a at H = 0 mT is a quantum effect and is not due to noise or heating in 
the measurement setup pnusH] . 
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Discussion 

In the literature there are no measurements of SCDs on junctions with a dominant second harmonic component 
in the CPR. Numerical simulations of the phase dynamics as a function of the damping parameter Q, the g factor 
and the temperature T give the conditions for which a double well potential effectively behaves in the escape process 
as a single well (for details see Supplementary Fig. 1, Supplementary Fig. 2 and Supplementary Note 1). Namely, 
for values of (5 — 10 SCDs with a single peak have been obtained for g < 2 or in the limiting case of pure second 
harmonic CPR (/i = 0). In this case the washboard potential changes its periodicity but assumes the form of a 
single well potential. Instead, for g > 2, two critical currents should be observed. In fact the heights of the two 
barriers approach each other when increasing the g factor, and the phase may be retrapped in both the potential wells 
with a finite probability, resulting in a bimodal switching distribution when counting many escape events [27ll29j (see 
Supplementary Fig. 3 and Supplementary Discussion). Since this is not observed, a pure second harmonic is the only 
possible explanation consistent with both measurements of magnetic field pattern[T^ (see Fig. Ic) and SCDs. 

Ico can be obtained by fitting the probability density of switching P{I) (red lines in Fig. 3b) in the thermal 
and quantum regime, Ico = 30.41 ± 0.05 p,A. A quite accurate value of the capacitance C can be obtained from the 
crossover temperature Tcross, which depends on Q, C and /co, see equation (1). By inserting the values of Ico and 
Q in the expression for Tcross we get C= 4.5 ± 0.9 pF. These values lead to ojp « 14 GHz. Nevertheless, as shown in 
Fig. 3c, the function a{g) in the equations for AU and ojp is a slowly varying function for g » 1, thus the junction 
parameters weakly depend on the g factor for high values of g. 

We expect in future additional insights coming from a comparative analysis with samples with lower SFE. For 
junction cross sections of about 50 /rm^ the Ic values for such junctions, are too high to be in the conditions to observe 
unambiguously the transition from the thermal to the quantum regime as commonly occurring also in standard S- 
Insulator-S iunctionspOl [55H55] . Only further advances in fabrication able to insert S-FI-S junctions in cavities 
and qubit architectures will give more refined feedback on the modes of the dissipative domains of the junction, and 
on the triplet component. However we can infer that the Q values of spin filter junctions are relatively higher than 
one would naively expect on the basis of the properties of the low SFE samples, which are characterized by higher 
values of Ic- Concerning the Q values of spin filter JJs, the reduction in Ic is compensated by the increase of Rn, as 
reported in Table 1. 

In conclusion, we have demonstrated the occurrence of MQT in NbN-GdN-NbN spin filter JJs. Spin filtering drives 
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the S-F-S junction in the underdamped regime and in the appropriate window of junction parameters to observe 
MQT. The switching current distributions, together with the period of Ic{H) modulation provides direct evidence for 
a pure second harmonic current-phase relation in the junction where MQT was observed. This is clear evidence of 
unconventional superconductivity, and it is possible that transport occurs by means of a pair of spin aligned triplet 
Cooper Dairs|16) which may suppress magnetic sources of decoherence |40L ITT] . Demonstration of macroscopic quantum 
phenomena in spin filter devices gives promise for their application in quantum hybrid circuits [TS], and also possibly 
as quiet memories. 


Methods 


Determination of spin filter efficiency. 

Spin filter efficiency (SFE) at a particular temperature is calculated from the RysT curve. SFE at any temperature 
below the Curie temperature {Tcurie) of GdN and above the superconducting transition temperature (Tc ) of NbN is 
calculated by defining SFE = where erf, crl are current densities of the up and down spin channels respectively. 

By invoking WKB approximation, we can derive SFE = tanh(cosh“^ {R*/R)) where R* is the measured value of 
resistance, and R is the value of resistance in the absence of spin filtering. R is estimated by fitting an exponential to 
the RvsT curve to temperatures above 50 K, and extrapolating the exponential to temperatures below the Tcurie of 
GdN. For calculating SFE below Tc of NbN, an exponential is fitted to the temperature dependence of SFE between 
TCurie and Tc and the fitted exponential is extrapolated to lower temperatures to obtain an estimate of SFE. 

Setup for switching current distributions. 

The SCDs have been measured by thermally anchoring the samples to the mixing chamber of He^/He"^ Oxford 
dilution refrigerator. The bias current is ramped at a constant sweep rate AI/At of about 2 mA s“^, the voltage is 
measured using a low noise differential amplifier and is fed into a threshold detector which is set to generate a pulse 
signal when the junction switches from the superconducting state to the finite voltage state. This signal is used to 
trigger a fast volt meter to record the value of the switching current. This procedure is repeated at least 10^ times at 
each temperature, which allows to compile a histogram of the switching currents. Filtering is guaranteed by a room 
temperature electromagnetic interference filter stage followed by low pass RG filters with a cut-off frequency of 1.6 
MHz anchored at 1.5 K, and by a combination of copper powder and twisted pair filters thermally anchored at the 
mixing chamber of the dilution refrigerator. 
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TABLE I: Device parameters. Spin filter efficiency has been determined at 4.2 K (see Methods for details), /co has been 
estimated from the fits of the switching current distributions and R„ from the IV curves. 

Barrier thickness (nm) Spin filter efficiency IcO (m^) (meV) IcRn (mV) 


1.5 

<30% 

820 

1700 

1.0 

1.7 

<30% 

280 

570 

0.9 

1.8 

60% 

120 

250 

0.7 

3.0 

90% 

30 

60 

0.1 
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FIG. 1: Spin filter Josephson junctions, (a) Thermal behavior of the resistance is a hallmark of spin filter NbN-GdN- 
NbN Josephson junctions (JJs). Semiconducting behavior has been observed at high temperatures, while below the Curie 
temperature of GdN (~ 33 K) the resistance decreases when decreasing the temperature. In non spin filter NbN-GdN-NbN JJs 
a semiconducting behavior has been observed down to the transition temperature of NbN. In panel b a picture of the device 
area is shown. The black arrow is the scale bar, 10 /rm. (c) Magnetic field pattern of a spin filter JJ in the range [-5 mT, 
5 mT], measured at 4 K. The blue curve is measured after nominal zero field cooling, the black and red curves are the magnetic 
pattern in the downward and upward direction of the magnetic field sweep, respectively, as indicated by the coloured dashed 
arrows. Blue squares indicate the values of the magnetic field, for which we present the measurements of switching current 
distributions, (d) Current-voltage characteristics of the same junction are shown as a function of the temperature, along with 
the voltage threshold (gray line) used for the switching current measurements. 
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FIG. 2: Phase dynamics with strong second harmonic component. Washboard potential in the case of a strong second 
harmonic component in the cnrrent-phase relation, for different values of the bias current 7 . Both panels a and b refer to the 
case g = 2, for which two potential wells are present. A and B indicate the lower and the higher potential well, respectively. 
The blue rectangle in a is magnihed in panel b, in which thermal activation (black arrow) and quantum tunneling processes 
(red arrow) from the well B at 7 = 0.9 are shown, ujp ( 7 , g) and A(7 ( 7 , g) are the plasma frequency and the height of the 


potential barrier respectively. 
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FIG. 3: Switching current distributions, (a) Measurements of switching current distributions (SCDs) from 4K down 
to 0.3 K. When increasing the temperature, the histograms move to lower values of switching currents Isw and the standard 
deviation a increases. Each histogram collects 10^ switching events, (b) Thermal behavior of the SCDs from 0.9 K down to 
20 mK. The red lines are fits of the probability density of switching, according to equation (2) in the case of a dominant second 
harmonic component in the current-phase relation (CPR). Below 100 mK the histograms overlap, indicating the transition to 
the quantum regime, definitely confirmed by the their dependence in magnetic field. In panel c the function a{g) is plotted 
(black curve) and the red line is the value of the function in the limiting case of a pure second harmonic (CPR). 
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FIG. 4: MQT in spin filter S-F-S Josephson junctions, (a) Thermal behavior of the mean switching current I mean (green 
up triangles, left axis) and of the standard deviation a (violet circles, right axis) in absence of magnetic field. The red dashed 
line indicates the crossover temperature between macroscopic quantum tunneling (MQT) regime (dark yellow background) 
and thermal activation (pink background). In the MQT regime both I mean and a saturate, (b) Comparison between thermal 
dependences of a at OmT (violet circles) and at 1.1 mT (orange right triangles). In presence of magnetic field, lower values 
of a have been measured and Teross is reduced, according to equation (1). Teross has been determined by the intersection of 
the curve in the thermal activation regime (dashed green lines) and the mean value of a in the MQT regime (black full 

lines). In (a) and in (b) the error bars have been determined by using the theory for propagation of statistical errors, from the 
equations defining the mean and the standard deviation of the switching histograms. 











